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ABSTRACT 


An  experimental  program  has  been  carried  out  to  determine  the  modes  of  raindrop 
breakup  in  air  after  passage  of  a  strong  shock  wave.  Relevant  dimensionless 
parameters  have  been  determined,  and  the  experimental  data  have  been  correlated 
in  terms  of  these  parameters.  The  resulting  correlations  have  been  used  to 
determine  raindrop  breakup  in  the  shock  layer  of  a  reentering  vehicle  passing 
through  a  rainstorm. 
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1.0  INTRODUCTION 


If  a  reentry  vehicle  flying  through  the  earth's  atmosphere  traverses  a  rainstorm, 
the  raindrops  striking  the  body  surface  at  high  speed  may  cause  serious  damage. 

This  damage  will  be  partially  alleviated,  however,  by  the  presence  of  the  flow 
field  around  the  body. 

Because  the  relative  velocity  between  the  rain  drops  and  the  flow  near  the  body 
is  high,  the  drops  are  decelerated  and  deflected  away  from  the  surface,  and  there 
is  a  tendency  for  the  drops  to  shatter  or  to  be  stripped  of  part  of  their  mass  as 
they  traverse  the  severe  environment  of  the  shock  layer.  In  the  latter  case,  the 
drops  may  be  broken  into  droplets  so  small  that  they  conform  to  the  local  stream¬ 
lines  in  the  shock  layer  and  thus  never  hit  the  body  surface. 

Seen  in  drop-fixed  coordinates,  the  stationary  drop  is  overtaken  by  the  bow  shock 
wave  of  the  vehicle,  and  the  flow  behind  that  shock  wave  (in  the  vehicle  shock 
layer)  fragments  the  drop.  This  view  of  the  problem  suggested  that  an  experimental 
study  of  the  shattering  could  be  conducted  in  a  shock  tube.  Such  experiments  are 
reported  herein  along  with  supporting  analysis  and  application  to  a  flight  case. 

Briefly,  the  purposes  of  the  study  were  as  follows:  ^irst,  to  determine  experi¬ 
mentally  the  time  after  passage  at  which  the  drop  was  ntirely  shattered;  second, 

to  determine  experimentally  the  drop  mass  as  a  function  >f  time  after  shock 
passage  and  until  complete  breakup;  third,  to  conduct  sufficient  analysis  to 
allow  general  application  of  the  experimental  data;  and  finally,  to  show  how  the 
experimental  and  analytical  results  could  be  applied  to  actual  flight. 

To  make  the  experiments  realistic,  the  tests  were  conducted  witn  shock  Mach 
numbers  of  3,  6,  9,  and  11  at  pressure  altitudes  of  sea  level,  20  kft,  and 
40  kft,  respectively.  The  corresponding  initial  pressures  were  760,  350,  and 
140  mm  Hg.  Water  drop  diameters  employed  varied  from  0.5  to  2.5  mm,  a  range 
that  is  representative  of  naturally  occurring  rain  drops.  The  experimental  data 
are  tabulated  in  the  appendix. 
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2.0  EXPERIMENTS 


2.1  EXPERIMENTAL  APPARATUS 

The  test  facility  employed  in  this  water  droplet  breakup  study  was  a  l^-inch- 
diamecer  shock  tube  capable  of  attaining  shock  Mach  numbers  up  to  12.  A  con¬ 
trolled  stream  of  water  droplets  was  introduced  into  the  tube  through  the  test 
section  located  at  the  end  of  the  driven  tube.  The  test  section  was  equipped 
with  viewing  ports  through  which  single-frame,  flash  X-ray  or  spark  shadowgraph 
photographs  were  obtained  at  times  varying  from  droplet/shock  wave  intercept  to 
disintegration.  A  schematic  of  the  facility,  shown  in  Figure  la,  indicates  the 
location  of  the  various  speed  and  pressure  monitoring  components.  A  closeup 
photograph  of  the  test  section  assembly  and  drop  generator  is  presented  in 
Figure  lb. 


|* - 30  in 

•  9-47M 


120  In 


Figure  la  SCHEMATIC  OF  TEST  FACILITY 


The  shock  tube  utilises  a  30- inch-long  driver  coupled  to  a  10-foot-long  driven 
section  between  which  is  located  a  steel-petal-type  diaphragm,  scored  to  rupture 
at  a  predetermined  pressure.  The  tube  is  operated  using  either  a  cold  hydrogen 
or  a  combustion  drive  to  cover  the  range  of  shock  Mach  numbers  required.  In  the 
cold  hydrogen  mode,  the  gas  was  fed  manually  at  a  very  slow  rate  into  the  driver 
until  the  diaphragm  ruptured.  Shock  speeds  up  to  Mach  3  were  attained  using 
this  mode.  To  attain  ahock  Mach  numbers  between  6  and  12,  it  was  necessary  to 
employ  the  combustion  drive  mods.  A  stoichiometric  mixture  of  hydrogen  and 
oxygen  with  a  70  percent  helium  dilution  was  ignited  by  an  electrically  exploded 
2-mil-dlameter  aluminum  wire  producing  driver  pressures  up  to  14,000  psia.  A 
driver  pressure  record  obtained  using  a  Klstler  piexo-elactrlc  transducer  is 
shown  in  Figure  2a  Indicating  a  slow  burning  rats.  Occasionally  very  fast  burning 
rates  were  obtained  as  shown  In  Figure  2b,  resulting  in  a  substantial  change  in 
shock  speed.  The  variation  in  burning  rates  was  attributed  to  the  random  mode 
in  which  the  wire  exploded  as  wall  as  the  degree  of  mixing  of  the  three  gaseous 
components.  Generally,  the  fast  burning  mode  resulted  in  unsatisfactory  test 
conditions. 

The  test  section  located  at  the  end  of  the  driven  tube  was  equipped  with  orthog¬ 
onal  ports;  one  set  for  Introducing  and  catching  the  droplet  stream,  the  other 
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Figure  2a  DRIVER  PRESSURE  RECORD  (SLOW  COMBUSTION) 
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Figure  2I>  DRIVER  PRESSURE  RECORD  (FAST  COMBUSTION) 


for  photographing  the  event.  Appropriate  seals  were  incorporated  into  these 
ports  to  maintain  the  subatmospheric  prerun  pressures  as  well  as  tc  contain  the 
high  transient  pressures  resulting  when  the  shock  wave  passes.  The  windows  are 
mounted  flush  with  the  inside  wall  of  the  test  section  to  minimize  distortion  of 
the  shock  wave.  The  water  inlet  and  outlet  ports  are  maintained  at  a  minimal 
diameter  for  the  same  reason.  Plexiglass  windows  were  used  as  viewing  ports  and 
were  replaced  frequently  since  the  surface  was  severely  deteriorated  by  the  high- 
temperature,  high-velocity  gas.  It  was  found  that  plexiglass  had  a  longer  life¬ 
time  under  these  conditions  than  glass  and  also  had  the  advantage  of  being  con¬ 
siderably  less  expensive  than  glass.  A  pressure  transducer  was  also  located  in 
the  test  section  to  monitor  the  pressure  rise  across  the  shock  and  consequently 
provide  a  double  check  on  the  shock  velocity.  A  typical  pressure  record  is  shown 
in  Figure  3.  The  slow  pressure  rise  indicated  is  the  result  of  the  slow  response 
of  the  pressure  transducer  to  the  step-function  pressure  increase. 

A  large  expansion  tank  is  located  at  the  end  of  the  tube  to  maintain  subatmos- 
Dheric  pressure  in  the  tube  after  the  test.  A  thin  Mylar  diaphragm  between  the 
tube  and  the  tank  enables  one  to  pump  the  tank  to  a  -much  lower  pressure  than  the 
tube.  The  primary  reason  for  using  the  expansion  tank  was  to  protect  the  droplet 
generator  from  long-duration  overpressure. 

The  electronics  used  to  synchronize  and  monitor  the  shock  location  and  spark 
light  source  or  flash  X-ray  triggering  is  sb jwn  schematically  in  Figure  4.  Four 
gages  located  in  the  wall  of  the  driver  tube  (see  Figure  1)  were  used  to  record 
the  speed  of  the  shock  waye.  The  first  gage  was  used  to  trigger  an  oscilloscope 
having  a  rastered  sweep.  Its  output  as  well  as  the  outputs  of  the  subsequent 
gages  were  displayed  on  the  scope  trace.  The  trigger  output  of  the  last  speed 
gage,  located  7  inches  from  the  droplet  stream,  was  used  to  trigger  the  spark 
gap.  A  delay  generator  was  used  lo  control  the  time  after  intercept  at  which 
the  trigger  occurred.  The  output  of  the  spark  was  also  fed  into  the  oscilloscope 
in  order  to  monitor  the  precise  location  of  the  shock  when  the  photograph  was  taken. 
A  typical  speed  record  is  shown  in  Figure  5a  and  the  corresponding  velocity  profile 
and  shock  location  graph  In  Figure  Sb..  In  this  shot  the  test  speed  was 
12,420  ft/sec. 

Two  types  of  gages  were  used  to  monitor  the  shock  velocity:  thin-film  heat  gages 
for  the  Mach  3  tests  and  ionization  gages  for  the  Mach  6  to  12  tests.  The  heat 
gages  consist  of  \  narrow  strip  of  platinum  painted  on  the  surface  of  a  Pyrex 
disk  which  was  mounted  flush  with  the  wall  of  the  driven  tube.  The  change  in 
temperature  resulting  from  the  passage  of  the  shock  was  sensed  and  recorded. 

The  fragility  of  this  gage  renders  it  impractical  for  the  high-speed  test;  conse¬ 
quently,  ionization  gages  were  used  for  the  Mach  6  to  12  tests.  This  gage  is 
essentially  a  modi! led  14  mm  spark  plug  mounted  flush  with  the  driven  tube  wall. 

The  gage  acts  as  the  open  end  of  a  make-circuit  which  is  closed  by  the  arrival 
of  the  ionized  shock  layer.  The  ionisation  gage  is  very  rugged  and  needs  oAly 
an  occasional  cleaning;  however,  it  is  too  insensitive  to  be  used  for  the  low 
Mach  number  tests. 

The  droplet  generator  system  was  based  on  a  technique  used  by  Dabora^,  in  which 
water  is  forced  into  a  cavity  (void  of  air  bubbles)  and  out  through  a  capillary 
as  shown  in  Figure  6.  (A  photograph  of  the  drop  generator  appears  in  Fig  e  lb.) 
Unlike  the  Dabora  System,  however,  this  system  was  designed  to  operate  in  vacuum 
tnd  under  very  high  transient  pressure  loading  (up  to  700  pals).  4  heavy-duty 
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Figure  3  TEST  SECTION  PRESSURE  RECORD 


Figure  4  SCHEMATIC  OF  ELECTRONICS 
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Figure  6  SCHEMATIC  OF  DROP  GENERATOR 


acoustic  speaker  driven  by  an  amplified,  variable  frequency,  sinusoidal  signal 
was  used  to  vibrate  a  3-mil-thick  brass  diaphragm  located  at  the  top  of  the  water 
chamber.  The  linkage  between  the  speaker  and  the  brass  diaphragm  was  attained 
by  cementing  a  Micarta  rod  to  the  speaker  diaphragm  and  coupling  it  to  a  copper 
rod,  soft-soldered  to  the  brass  diaphragm.  An  air-bleed  hole  was  drilled  through 
the  copper  rod  to  allow  the  purging  of  the  chamber  of  air  after  each  test.  The 
capillary  was  epoxy  adhesive-bonded  into  the  end  of  the  water  chamber  opposite 
the  diaphragm.  The  cavity  assembly  was  clamped  together  with  three  1/4-20  bolts 
which  compressed  the  two  0-ring  seals  at  each  end  of  the  droplet  generator. 
Distilled  water  was  driven  into  the  chamber  using  nitrogen  gas  pressure  as  shown 
in  Figure  6.  The  water  flow  into  the  chamber  was  controlled  by  a  manually  ad¬ 
justed  micrometer  valve.  The  flow  rate  used  for  all  capillary  sizes  was  that 
rate  which  was  adequate  to  maintain  a  chamber  back  pressure  sufficient  to  estab¬ 
lish  a  laminar  water  jet  at  the  orifice  of  the  capillary.  Three  capillary  sizes 
were  used  (0.007-,  0.018-,  and  0.060-inch  I.D.)  for  three  drop  nominal  diameters, 
0.5,  1.0,  and  2.1  mm,  respectively.  The  critical  frequency  at  which  the  brass 
diaphragm  must  be  vibrated  to  sustain  a  stream  of  droplets  hf-ving  uniform  size, 
spacing,  and  shape  was  determined  experimentally  for  each  capillary.  The  droplet 
stream  was  collected  in  a  flask  attached  to  the  underside  of  the  test  section 
which  was  drained  after  each  test. 

Two  optical  systems  were  used  during  the  course  of  the  tests,  namely,  a  single¬ 
frame  spark  shadowgraph  and  a  single-frame  flash  X-radiograph.  The  former 
optical  system  is  shown  schematically  in  Figure  7a.  A  spark  light  source  having 
an  effective  duration  of  0.3  Msec  is  used  to  backlight  the  event.  The  light 
emitted  by  the  spark  gap  is  collimated  by  an  f/2.5  aero-ektar  lens  prior  to  pass¬ 
ing  through  the  object  plane.  The  parallel  beam  of  light  is  subsequently  passed 
through  the  objective  lens,  which  serves  a  twofold  function.  First,  it  focusses 
the  droplets  on  the  film  plane,  and  second,  it  images  the  spark  gap  on  the  camera 
iris  which  is  closed  down  to  its  smallest  opening.  This  restricting  aperture 
limits  the  amount  of  stray  radiation  (primarily  from  the  luminous  air  behind  the 
shock)  that  can  strike  and  thereby  fog  the  film.  A  further  reduction  in  the 
fog  level  was  achieved  by  introducing  an  interference  filter  (4220  A)  between 
the  camera  and  the  objective  lens.  The  bandpass  of  the  filter  was  selected  to 
pass  the  radiation  from  the  spark  while  rejecting  the  radiation  from  the  heated 
air. 

A  Fexitron  flash  X-ray  system  utilizing  a  model  524  X-ray  tube  operating  at  a 
voltage  of  105  kv  was  used  in  the  X-radiograph  system  as  shown  in  Figure  7b. 

The  effective  source  size  is  less  than  1  mm  and  has  a  duration  of  30  nsec.  A 
diverging  beam  of  X-rays  is  emitted  by  the  X-ray  source  which  passes  through  the 
test  section  onto  the  film,  located  as  close  to  the  object  plane  as  possible. 
Polaroid  X-ray  film,  type  TLX,  was  used  because  of  its  fine  grain  and  high 
sensitivity  (speed).  The  magnification  of  the  object  is  determined  by  the  ratio 
of  the  source-object  distance  to  the  source-film  distance. 

Resolution  of  2  ram-diameter  water  droplets  with  the  above  system  could  not  be 
attained  without  the  addition  of  a  lead  salt  to  the  water.  Adequate  increase  in 
X-ray  absorption  was  obtained  by  dissolving  lead  acetate  trihydrate  to  the  point 
of  saturation.  This  solution  was  carefully  filtered  to  remove  undissolved 
material  and  stored  in  a  sealed  container  to  prevent  evaporation  of  liquid  and 
subsequent  precipitation  of  solute.  The  effect  of  doping  the  water  with  lead 
acetate  was  to  increase  the  density,  viscosity,  and  surface  tension  of  the 
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solution  above  those  for  pure  water.  Specifically  the  ratios  of  the  density, 
viscosity,  and  surface  tension  of  the  solution  to  pure  water  values  are  1.44, 

1.82,  and  1.05  respectively.  As  will  be  seen  later,  proper  non-dimensionaliza- 
tion  of  the  experimental  results  renders  the  data  obtained  with  the  doped  water 
directly  comparable  to  data  obtained  using  pure  water. 

A  schematic  view  of  the  drop  shattering  experiment  in  room  fixed  coordinates  is 
shown  in  Figure  8.  A  drop  of  initial  diameter  and  density,  D0  and  pi  respectively, 
is  shown  prior  to  the  passage  of  a  shock  moving  with  a  speed  us  into  an  initial 
pressure  p.  Gas  conditions  behind  the  shock  are  indicated  by  subscript  2.  Time 
is  measured  starting  from  the  instant  the  shock  crosses  the  drop.  The  distance 
x  is  the  coordinate  of  the  drop  relative  to  its  initial  location,  and  unsub- 
scripted  u  is  the  drop  speed,  dx/dt. 
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Figure  8  SCHEMATIC  OF  EXPERIMENT 
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2.2  SPARK  SHADOWGRAPH  DATA 

Spark  shadowgraph  data  vara  obtained  for  three  reasons s  First,  we  wished  to 
determine  accurately  the  distortion  history  of  the  water  drops  at  timas  immediately 
after  shock  passage.  These  data  were  necessary  to  substantiate  the  first  step  in 
the  theoretical  analysis  of  drop  behavior.  Second,  it  was  necessary  to  show  that 
the  doped  drops  behaved  the  same  as  the  undoped  drops  in  properly  non-dimension- 
all zed  time.  Finally,  the  spark  shadowgraphs  provided  a  simpler  method  than  X  rays 
to  obtain  an  approximate  idea  of  when  the  drops  were  shattered  under  the  most 
severe  test  condition,  M,  -  11,  p  ■  141  torr.  Studies  already  reported  indicated 
approximately  what  the  time  to  breakup  was  under  low  Mach  number  test  con¬ 
ditions.  2, 3, 4 
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Immediately  after  passage  of  the  shock  wave  over  the  drop,  the  flow  over  the 
drop  sets  up  a  pressure  differential  between  the  drop's  stagnation  point  (or  pole) 
and  its  shoulder  (or  equator).  In  response  to  this  pressure  difference  the  drop 
expands  normal  to  the  flow  direction  and  contracts  in  the  direction  parallel  to 
the  flow  direction.  This  shape  change  is  illustrated  in  Figures  9a  through  13a. 
These  shadowgraphs  were  obtained  for  distilled  water  drops  with  a  1.03  mm  initial 
diameter  overtaken  by  a  Mach  6  shock  at  an  initial  pressure  of  350  torr.  Fig¬ 
ure  9a  shows  a  drop  prior  to  shock  arrival,  while  the  subsequent  shadowgraphs 
correspond  respectively  to  times  after  shock  passage  of  approximately  4,  9,  15, 
and  25  ^sec. 

Ranger  and  Nicholls^  have  reported  that  for  their  data  at  shock  Mach  numbers  be¬ 
low  4,  the  drop  distortion  r/r0  (where  r  is  the  radius  of  the  drop  at  its  equator) 
could  be  correlated  against  a  dimensionless  time  after  shock  passage  given  by 
T  =  tu  2  b  /  Do  where  b  =  p  2  /  pL  ■ 

To  obtain  accurate,  general  early-time  distortion  histories,  data  were  acquired 
in  the  present  tests  over  the  entire  experimental  range;  i.e.,  for  Ms  from  3  to 
11,  D0  from  0,5  to  2.5  mm,  p  from  760  to  141  torr,  and  doped  or  undoped.  The 
results  are  shown  in  Figure  14,  where  t/t0  is  plotted  against  dimensionless 
time  T. 

Since  the  time  scale,  d0 /ujb H  ,  varies  by  a  factor  of  10  in  these  experiments, 
Figure  14  indicates  that  T  is  the  appropriate  dimensionless  time,  and  also  in¬ 
dicates  that  r(T)/ r0  is  concave  upward  at  early  time.  The  significance  of  these 
data  with  regard  to  the  analysis  will  be  discussed  later. 

Similar  results  for  one  test  condition,  but  for  later  time  are  shown  in  Figure  15. 
Here  distortion  data  for  drops  saturated  with  lead  acetate  are  compared,  in  dimen¬ 
sionless  time,  against  results  obtained  by  Reinecke  and  McKay ^  at  the  same  test 
condition  (M=6),  Figure  15,  along  with  the  two  doped  points  in  Figure  14,  in¬ 
dicates  that  the  density  difference  between  doped  and  undoped  drops  is  adequately 
accounted  for  by  the  appearance  of  p^/i  in  T. 

The  spark  shadowgraph  system  was  also  employed  to  determine  the  approximate  time 
required  to  shatter  the  drops  at  the  most  severe  test  conditions,  Ms  “  11,  p  ■  141. 
These  tests  indicate  a  dimensionless  time  to  breakup,  Tt>»  °f  about  3,  distinctly 
below  the  values  of  5  to  6  reported  at  Mach  2  to  3.5,  also  on  the  basis  of  shadow¬ 
graph  data,  by  Ranger  and  Nicholls.4 

Finally,  Figure  11a  indicates  a  spark  shadowgraph  of  a  drop  at  Mach  6,  30  psec 
after  shock  passage  (T"  1.28).  By  this  time  a  substantial  amount  of  water  has 
been  stripped  from  the  equator  of  the  drop  and  has  entered  itB  wake,  rendering 
the  wake  opaque  to  visible  light.  Comparison  of  Figure  11a  with  the  undisturbed 
drop  shown  in  Figure  9a  indicates  that  simple  mass  continuity  precludes  the  pos¬ 
sibility  that  the  entire  shadow  in  Figure  11a  contains  only  water.  Clearly  the 
remaining  unstripped  mass  of  the  original  drop  is  confined  to  a  lens-shaped  volume 
at  the  front  of  the  shadow  in  Figure  11a.  How  much  water  remains,  however,  cannot 
be  determined  from  the  ahadowgraph.  Instead,  the  flash  X-ray  system  was  employed 
to  measure  the  remaining  unatripped  mass  of  the  drop. 
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UNDISTURBED  DROP  PRIOR  TO  SHOCK  ARRIVAL 
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b.  ISOPHOTE  OF  RADIOGRAM 


Figure  14  EARLY-TIME  DROP  DISTORTION 
HISTORY 


Figure  15  A  COMPARISON  OF  DROP  DISTORTION 
HISTORIES  FOR  DOPED  AND 
UNDOPED  DROPS 


2.3  X-RADIOGRAM  DATA 

Where  a  pulse  of  X-rays  traverses  a  homogeneous,  partially  X-ray-transparent 
body,  the  ratio  of  transmitted-to-incident  X  rays  or  photons  is  proportional 
—  Ci  z 

toe  ,  where  *  is  the  local  thickness  of  the  body  along  the  X-ray  path.  More¬ 
over,  in  the  "linear"  response  region  of  an  X-ray  plate,  the  photographic  density 
is  proportional  to  the  logarithm  of  the  flux  or  number  of  photons  striking  the 
film.  Combining  these  relations  immediately  shows  that  the  density  of  the  X-ray 
plate  is  given  by  C1+C2*,  where  Cj  is  simply  the  background  density  of  the  film, 
i.e.,  where  ■  0.  Thus,  if  the  path  length  s  is  known  at  any  point  corresponding 
to  the  X-ray  image,  C2  can  be  calculated  and  the  path  length  elsewhere  can  be 
computed  from  the  local  relative  film  density.  Even  if  the  film  response  is  not 
strictly  linear,  the  relation  between  film  density  and  path  length  can  be  ob¬ 
tained  experimentally  from  the  X  radiogram  of  a  body  of  known,  varying  thickness. 

Now  consider  a  drop  of  water  containing  a  lead  salt  in  solution.  X  rays  traversing 
the  drop  will  be  absorbed  by  the  lead  ions.  Assuming  that  the  X  rays  are  uniformly 
dispersed  in  the  water,  the  film  density  in  the  resulting  radiogram  will  be  propor¬ 
tional  to  the  amount  of  water  along  the  X-ray  path.  Then  a  scanning  microdensit  >- 
meter  may  be  used  to  measure  the  distribution  of  film  density  over  the  drop  image 
in  the  X-ray  plate  by  plotting  lines  of  constant  film  density  and  the  drop  mass 
or  volume  obtained  from  these  "lsophote"  diagrams  by  integration,  in  tha  same 
manner  that  the  volume  of  a  mountain  nay  be  calculated  from  a  contour  map. 
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This  method  was  employed  to  measure  the  mass  remaining  in  liquid  drops  as  a 
function  of  time  after  shock  passage.  The  scanning  microdensitometer  isophote 
of  a  doped,  spherical,  2.12-mm-diameter ,  undisturbed  water  drop  is  shown  in 
Figure  9b,  The  slight  apparent  oblateness  of  the  drop  is  due,  not  to  real  non¬ 
sphericity,  but  to  a  slightly  different  magnifications  in  the  vertical  and 
horizontal  scan  of  the  microdensitometer.  The  increasing  film  density  and  cor¬ 
responding  sphere  secant  is  indicated  in  the  isophote  by  the  digitized,  nine- 
unit.  color  code  of  dotted  red,  solid  red,  white,  dotted  blue,  solid  blue,  white, 
dotted  green,  solid  green,  and  white.  The  film  density  versus  water  path  length 
calibration  obtained  from  this,  and  a  similar  isophote,  is  shown  in  Figure  16. 

The  calibration  points  were  obtained  by  scanning  along  orthogonal  diameters  on 
each  of  two  different  radiograph  of  different  drops.  Thus  the  data  represent 
four  scans.  Although  it  is  not  apparent  on  the  distorted  scale  of  Figure  16, 
the  calibration  indicated  a  nearly  linear  X-ray  plate  response  except  near  the 
image  edges  where  the  finite  size  of  the  microdensitometer  aperture  distorted 
the  calibration.  To  check  the  repeatability  of  film  responses  from  plate  to 
plate,  of  the  drop  diameter,  of  the  film  and  densitometer  magnification,  etc., 
the  calibration  was  used  to  calculate  the  volume  of  an  undisturbed  drop  X-rayed 
late  in  the  test  program.  The  mass  of  the  drop  was  computed  to  be  0.93  of  the 
calibration  drop,  indicating  a  very  satisfactory  level  of-  overall  repeatability. 

The  integrations  to  obtain  drop  mass  were  carried  out  on  ah  IBM  360  computer 
using  as  input  a  tape  of  the  microdensitometer  output  corresponding  to  the 
isophotes  and  digitized  over  64  shades  of  intensity. 

Typical  isophotes  of  doped,  X-rayed  drops  at  various  times  after  shock  passage 
are  shown  in  Figures  10b  through  13b  and  correspond  in  dimensionless  time,  T,  to 
the  respective  shadowgraphs.  The  aerodynamics  parameters  were  the  same  in  these 
X-ray  tests  as  is  illustrated  in  Figures  10a  through  13a;  i.e.,  Mach  6  at  350  torr. 
However,  the  initial  drop  diameter  for  the  X-ray  tests  was  2.12  mm  (this  diameter 
yielded  significantly  better  X-ray  resolution  than  1  mm  drops),  and  the  water  was 
saturated  with  lead  acetate,  which  increased  the  density,  viscosity,  and  surface 
tension.  As  expected.  Figures  10,  11,  and  12  indicate  that  the  remaining  un¬ 
stripped  water  is  confined  to  a  lenticular  volume  at  the  front  of  the  shadow¬ 
graph  image  and  that  the  stripped  mass  is  very  tenuously  distributed.  Figure  13 
shows  the  drop  at  a  fairly  late  state  of  disintegration.  In  this  particular  test 
it  has  broken  into  two  smaller  drops.  The  shadowgraph  in  Figure  13a  was  chosen 
from  several  available,  because  it  was  geometrically  compatible  with  the  isophote. 
The  drop  separation  into  two  distinct  droplets  was  not  generally  seen  in  the  tests. 
In  fact,  as  will  be  discussed  later,  the  final  disintegration  process  results  in 
a  random  and  unpredictable  geometry. 

From  such  isophotes,  the  drop  mass  histories  shown  in  Figures  17,  18,  and  19  were 
obtained.  At  the  Mach  3  condition  (see  Figure  17),  the  drop  mass  decreases 
smoothly,  and  indicates  a  time  to  breakup  of  about  150  psec.  At  Mach  6  (see 
Figure  18) ,  the  drop  mass  decreases  smoothly  to  about  40  percent  of  its  initial 
value,  at  which  time  (about  50  microseconds)  the  mass  decreases  abruptly. 

In  the  Mach  11  case  (see  Figure  19)  the  drop  image  on  the  X-ray  plates  abruptly 
dropped  below  the  background  density  after  about  25  psec.  Several  attempts  were 
made  to  obtain  X-ray  images  at  about  25  psec.  These  tests  apparently  yielded 
either  dark  X-ray  plates  or  Images  whose  computed  mass  exceeded  70  percent  of 
the  original  mass.  Those  points  showing  either  70  percent  or  more  or  a  dark 
plate  may  be  seen  in  Figure  19.  This  result  is  apparently  caused  by  an  extremely 
abrupt,  catastrophic  disintegration  of  the  drop  at  about  25 psec  after  shock 
passage. 
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Figure  16  X-RAY  PLATE  CALIBRATION 


3.0  CORRELATION  AND  INTERPRETATION  OF  EXPERIMENTAL  DATA 


The  experimental  data  reported  in  Section  2.0  has  been  correlated  in  terms  of 
appropriate  nondimensional  parameters.  The  following  material  describes  the 
derivation  of  the  nondimensional  parameters  from  theoretical  considerations  and 
the  correlation  of  the  data. 

3.1  EXPERIMENTAL  DROP  TRAJECTORIES 

e 

Wc  will  here  simply  cite  results  already  reported  in  the  literature.  The  tra¬ 
jectory  of  a  drop  is  strongly  affected  by  the  deformation  of  the  drop,  which  is 
caused  by  the  suddenly  imposed  pressure  differential  between  the  equator  and  the 
poles.  This  deformation  changes  the  drag  coefficient  of  the  drop,  and  the  drop 
trajectory  becomes  difficult  to  compute.  The  stripping  of  mass  from  the  drop  is 
also  coupled  to  the  trajectory  computation.  Under  such  circumstances  it  becomes 
simpler  to  determine  the  drop  trajectory  experimentally. 

The  tiajectory  of  a  deforming  drop  may  be  correlated  in  terms  of  a  dimensionless 
time  based  on  the  time  constant  of  the  deformation.  The  magnitude  of  the  pressure 
differential  is  given  by  the  dynamic  pressure  1/2  -u2 2 actin8  on  t^e  drop. 

Dimensional  analysis  shows  that  the  characteristic  time  associated  with  deforming 
the  drop  by  an  amount  on  the  order  of  the  original  diameter,  l)Q  ,  in  some  direc¬ 
tion  is  given  by 


Hence  we  can  define  a  nondimensional! zed  time,T  ,  by  the  equation 


Let  X-  x  0o  be  the  distance  in  original  drop  diameters  that  the  drop  has  moved 
in  the  shock  tube  at  time  i  .  The  drop  trajectory  data  compiled  at  Avco  chn  be 
correlated  by  the  relationship^ 

X  *  O  R  TJ,  (3) 

which  corresponds  to  a  constant  acceleration  trajectory  with  an  average  drag 
coefficient  of  about  2.1. 

3.2  DROP  DEFORMATION 

The  complete  analysla  of  the  deformation  of  a  waterdrop  following  the  passage  of 
a  shock  wave  is  a  formidable  problem  far  beyond  the  scope  of  the  present  study. 
Instead  of  attempting  a  complete  numerical  or  analytical  solution  to  the  equations 
of  motion  for  tbu  drop,  we  have  explored  two  approximate  analyses  in  order  to 
gain  some  insight  into  the  general  character  of  the  deformation  and  the  governing 
parameters. 
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In  the  first  approach,  the  initial  deformation  of  the  drop  was  considered.  Pro¬ 
vided  attention  is  confined  to  the  period  immediately  following  the  shock  passage 
(the  possible  length  of  such  a  period  will  be  explored  subsequently),  the  boundary 
of  the  dror>  may  be  approximated  as  the  initial  sphere,  and  the  solution  to  the 
equations  of  motion  can  be  developed  in  spherical  harmonics.  Let  cf>  be  the  velo¬ 
city  potential  for  the  motion,  satisfying  the  Laplace  equation  in  spherical 
coordinates 

42*  .  „i„»  |  .  £  ("*  w)  '  °  '  (4) 

where  r  is  the  radial  coordinate  and  6  is  the  azimuth  angle  measured  from  the 
pole  of  the  drop  (i.e.,  the  axis  of  symmetry  of  the  flow). 

The  general  solution  to  Equation  (4)  is 


<f>  -  ^  A„  (t)  rn  P„  (0)  , 

n 

where  Pn  are  the  Legendre  polynomials  in  n  =  cos 

p0  =  l.Pj  fi,  p2  =  J  (3  ft2  -  1),  P3  =  J  (5M3  -3m) . 


(5) 


It  can  be  shown  that  an  impulsive  pressure  ~  is  related  to  an  imposed  value  of  <6 
on  the  boundary  of  the  drop  by  the  equation  <4  m  &J/pL  +  constant,  where  is  the 
water  density.  Since  in  our  case  the  "impulsive  pressure"  is  actually  applied 
continuously,  we  have  ~  =  pt  .  In  addition,  the  arbitrary  constant  just  referred 
to  can  be  taken  to  vanish  without  loss  of  generality.  We  then  have  at  the  sur¬ 
face  of  the  drop  r  *  r 

O 


*  -  5  V'K  pnw  =  •  (6) 

Pl 

which  is  sufficient  to  determine  the  coefficients  for  e  given  pressure  distri¬ 
bution. 

This  approach  has  been  taken  by  Burgers  (referenced  in  Engel2)  in  the  case  of  a 
symmetrical  pressure  distribution  of  the  form 

P  =  Pe  +  (ps  -  pe)  cos2  0  (75 
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such  as  would  be  obtained  in  subsonic  lnviscid  flow.  He  then  obtains 
1  PS  *  2Pe 


2  Ps  "  Pe  t 


A1  =  A3  =  a4  -  A5  =  .  .  .  =  0  . 


(8) 


The  initial  deformation  of  the  drop  can  be  obtained  from  this  solution  by  ex¬ 
pressing  the  radial  velocity  component  v  in  terms  of  the  velocity  potential: 


V  = 


dr 

dt 


dd> 

dr 


A„rn 


(9) 


By  integrating  in  time,  we  get  the  relative  deformation  of  the  contour  at  any 
point  on  the  periphery  of  the  drop: 


(10) 


In  terms  of  the  Burgers'  solution,  this  displacement  becomes  at  the  poles  of  the 
drop 


2  Ps  ~  Pe  t2 


and  at  the  equator  of  the  drop 


1  Ps-Pe  t2 


In  her  report,  Engel2  compares  Burgers'  result  with  her  experimental  data  on  the 
equatorial  dilatation  of  vaterdrops.  There  appears  to  be  a  numerical  error  in 
Engel's  computations  of  the  Burgers'  solution,  however,  and  the  epparent  agree¬ 
ment  of  experiment  with  theory  shown  in  her  report  is  spurious.  We  have  reduced 
Engel's  data  and  Burgers'  theory  in  terms  of  the  dimensionless  time 
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where  u2  is  the  initial  gas  velocity  relative  to  the  drop,  p2  is  t'ie  8as  density 
and  D0  is  the  initial  drop  diameter. 

Burgers'  solution  then  becomes  simply 


D  2 
8  s  —  -  1 

D„ 


-  T2 

3  P 


(14) 


for  the  relative  equatorial  dilatation,  where 


C 


P 


Ps  *  Pe 

1/2  p2  u| 


(15) 


The  solution  is  plotted  in  Figure  20  for  the  subsonic  incompressible  value 
C  =  9/4,  and  it  can  be  seen  that  although  the  experimental  data  are  fairly  well 
correlated  in  terms  of  the  dimensionless  time,  the  agreement  with  Burgers'  solu¬ 
tion  is  not  satisfactory. 

It  is  worth  noting  that  recently  Ho&  has  derived  much  more  complete  solutions  to 
the  drop  deformation  equations,  and  that  his  solutions  are  in  general  agreement 
with  Burgers'  simplified  analysis  for  the  initial  deformation  assuming  an  in- 
viscid  pressure  distribution. 


In  an  attempt  to  determine  the  effect  of  the  assumed  pressure  distribution  on  the 
analytical  result,  we  derived  the  solution  for  the  initial  deformation  for  a  case 
in  which  the  pressure  had  the  same  (subsonic)  variation  over  the  windward  side  of 
the  drop,  but  was  constant  over  the  leeward  side  of  the  drop  rather  than  being 
symmetrical  with  respect  to  the  equatorial  plane.  Such  a  pressure  distribution 
should  be  a  rough  approximation  to  the  actual  distribution  on  a  drop  with  a  wake. 
Fourier-analyzlng  this  distribution  and  truncating  at  the  fourth  term,  we 
obtained 


P 


Pe  + 


co*  8  * 


1 

—  co*  28  + 


8 

15* 


cos 


*) 


(16) 


to  within  a  percent  or  so  of  the  postulated  distribution.  Substituting  in  the 
above  analysis,  we  determined  that  the  corresponding  equation  for  the  dilatation 
of  the  drop  equator  is 


1  P*  ’  Pe  t2 


half  the  Burgers'  result,  and  even  more  out  of  line  with  the  experimental  data. 

Figure  21  shows  the  deformation  of  the  drop  prodicted  by  this  analysis  after  an 
arbitrary  small  time.  It  may  be  observed  that  the  drop  tends  to  flatten  con¬ 
siderably  more  on  the  windward  than  on  the  leeward  side,  and  that  the  point  of 
maximum  lateral  deformation  la  slightly  forward  of  the  original  equator.  The 
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DIMENSIONLESS  TIME,  T 


'Figure  20  DROP  DEFORMATION:  COMPARISON  OF  ENGEL'S  DATA  AND 
BURGERS'  SOLUTION 


aoisw 

Figure  21  DROP  FLATTENING  WITH  SIMULATED  WAKE  FLOW 
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maximum  lateral  dilatation  is  close  to  the  dilatation  at  the  original  equator, 
however,  and  the  observation  just  stated  still  holds  true. 

At  this  point,  it  was  concluded  that  the  analysis  of  the  initial  deformation  of 
the  drop  was  simply  not  applicable  to  the  observed  experimental  behavior.  The 
reason  is  probably  either  that  the  physical  time  scale  of  the  observed  deformation 
is  much  larger  than  the  time  scale  which  describes  the  region  of  validity  of  the 
analysis,  or  that  the  assumption  of  "impulsive  pressures"  is  somehow  invalid.  In 
either  case,  the  analysis  is  mathematically  correct  but  physically  irrelevant. 

Our  second  approach  to  a  solution  to  the  deformation  problem  was  more  successful. 
In  this  approach  we  discarded  the  possibility  of  an  exact  mathematical  solution, 
no  matter  how  restricted  the  region  of  validity,  and  considered  the  drop  to  be 
disk-shaped  rather  than  nearly  spherical.  The  reason  for  making  this  assumption 
was  that  a  simple  approximate  solution  could  then  be  obtained  for  all  times. 

Let  R  be  the  radial  coordinate  in  a  cylindrical  coordinate  system  with  axis  along 
the  axis  of  symmetry  of  the  drop,  and  let  v  be  the  corresponding  fluid  velocity 
component  in  the  drop.  Our  primary  assumption  is  that  the  flow  variables  do  not 
depend  on  the  axial  coordinate  but  are  (in  some  approximate  sense)  uniform  from 
front  to  back  of  the  drop  at  any  radius  R.  This  assumption  cannot  be  justified 
a  priori,  but  it  should  be  borne  in  mind  that  the  purpose  of  this  analysis  is  to 
get  a  rough  idea  of  the  deformation  of  the  drop  rather  than  a  detailed  solution. 
The  radial  equation  of  motion  in  the  drop  then  becomes 


av  dv  1  dp 

dt  dR  dR 


(18) 


We  now  integrate  this  equation  from  the  axis  R  ■  0  to  the  rim  R  ■  i  : 


d 

dt 


vdR  + 


■j  v2(r) 


Ps'Pe 


(19) 


The  outflow  velocity  la  assumed  to  vary  as  some  power  of  R : 


Substituting  In  Equation  (19)  and  Integrating,  we  obtain 
t  d*r  3  +  k  /  dr  \2  P«  “  Pe 

f7k  +  2<i+k)  \dT;  "  PL  • 

or  In  dimensionless  variables 
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3+  k 
2(1  +  k) 


(21) 


1  »  <3  d2  <5 


Equation  (21)  has  been  integrated  numerically  for  various  values  of  the  exponent 
k.  The  results  shown  on  Figure  22  demonstrate  that  the  solution  is  not  a  very 
sensitive  function  of k  ,  and  we  shall  use  the  solution  for  k  =  1  as  the  refer¬ 
ence  solution. 

The  main  characteristics  of  the  reference  solution  are  that  it  initially  follows 
the  "transient"  solution 

5  -  2  CP  T2  (22) 

(which  incidentally  is  three  times  larger  than  the  Burgers'  solution),  but  that 
it  very  rapidly  approaches  the  "steady-state"  solution 

S  -  yflC ~  T  +  const.  (23) 

Essentially  the  steady  state  is  reached  by  T  »  1/2.  We  may  infer  from  this  that 
the  initial  solutions  discussed  previously  are  probably  invalid  because  their 
region  of  validity  is  too  restricted. 

The  reference  solution  has  been  plotted  in  Figure  20,  and  it  is  seen  to  offer 
somewhat  better  agreement  with  Engel's  data  than  with  the  Burgers'  solution. 

Its  agreement  with  the  higher-speed  data  of  Ranger  and  Micholls^  and  the  present 
study,  however,  is  amazingly  good,  as  shown  in  Figure  23.  It  may  be  observed 
that  the  model  of  a  disk-shaped  drop,  while  rather  crude  for  a  nearly  spherical 
drop,  is  perhaps  not  unrealistic  at  later  times  when  the  drop  has  been  con¬ 
siderably  flattened. 

3.3  CATASTROPHIC  MODE 


The  dominant  effect  in  the  catastrophic  mode  is  seen  to  be  the  rapid  growth  of 
surface  waves  on  the  windward  face  of  the  drop.  Taylor7  pointed  out  that  when 
the  direction  of  acceleration  of  a  fluid  interface  is  from  the  less  dense  fluid 
to  the  more  dense,  waves  on  the  interface  will  be  unstable  and  will  grow  ex¬ 
ponentially.  Small  waves  or  irregularities  are  always  present  in  any  real 
physical  situation,  but  when  they  are  unstably  amplified,  they  may  grow  so 
large  that  they  tear  the  drop  apart. 

Taylor's  theory7  for  the  growth  of  the  waves  is  strictly  applicable  only  to 
small  disturbances,  but  it  may  be  interpreted  as  providing  the  relevant  param¬ 
eter  for  correlating  the  experimental  data  even  though  the  wave  amplitudes  be¬ 
come  very  large  (at  the  point  of  disintegration,  the  wave  amplitudes  are 
comparable  in  magnitude  to  the  diameter  of  the  drop  itself).  The  theory  predicts 
the  following  characteristic  amplification  time  for  the  wave  amplitudes: 
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(24) 
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The  corresponding  nondimensional  time  is 


I  4 


1  4  ,,1  2  nl  4 

pi  '-‘I  uo 


(25) 


Hence  the  theory  suggests  the  following  correlation  parameter  for  the  catastrophic 
mode: 


l>2V  2  Do 
Vi  - 


(26) 


which  is  the  Weber  number. 


We  now  turn  to  the  relevant  experimental  data.  The  variation  of  drop  mass  with 
tine  shown  in  Figures  17,  18,  and  19  corresponding  to  tests  at  Mach  3,  6,  and  11, 
respectively,  indicates  catastrophic  breakup  is  occurring  early  (with  respect  to 
stripping  effects)  in  the  Mach  11  tests,  at  intermediate  time  at  Mach  6,  and  late 
time  (if  at  all)  at  Mach  3.  In  view  of  the  analysis  given  above,  it  is  appro¬ 
priate  to  plot  these  breakup  times  versus  the  Weber  number,'*  .  The  times  to 
breakup  in  the  three  cases  are  150,  50,  and  22  msec,  respectively,  or  in  the 
dimensionless  form,  2.7,  2.0,  and  1.4.  These  values  of  Tb  are  plotted  versus 
«•  in  Figure  24,  and  are  seen  to  be  well  correlated  by  the  equation 

Tb  45  W 1  4  (27) 

Thus  the  theoretical  analysis  of  the  catastrophic  mode  appears  well  verified. 
Moreover,  it  appears  that  even  the  lowest  Mach  number  data  were  subject  to 
catastrophic  breakup,  but  only  after  a  substantial  amount  of  stripping  had 
occurred. 

Further  verification  of  this  dependence  of  Tb  on  Weber  number  is  available  in 
the  data  of  Reinecke  and  McKay. ->  They  report  high  Weber  number,  catastrophic 
breakup  data  on  the  basis  of  their  spark  shadowgraphic  measurements  of  time  to 
breakup.  These  four  relevant  data  are  also  plotted  in  Figure  24  and  are  reason¬ 
ably  well  correlated  by  W-l/\  Thus  the  spark  shadowgraphic  method  of  deter¬ 
mining  breakup  time  appears  to  give  reliable  qualitative  results,  but  generally 
to  overestimate  breakup  time  by  about  a  factor  of  two,  at  least  in  this  high 
Weber  number  regime. 

The  values  of  Weber  number  employed  in  these  calculations  were  based  on  the  room- 
temperature  value  of  the  water  surface  tension  .  It  may  be  objected  that  the 
surface  of  the  drop  is  raised  to  much  higher  temperatures  by  aerodynamic  heating 
in  the  shock-tube  experiments,  so  that  a  lower  value  of  a.  should  be  used.  We 
note,  however,  that  the  thermal  layer  in  the  water  drop  is  extremely  thin  under 
transient  experimental  conditions,  and  that  the  choice  of  an  appropriate  value 
for  depends  (among  other  things)  on  considerations  of  the  interaction  of  this 
thermal  layer  with  the  growing  surface  waves  themselves.  Such  considerations  are 
beyond  the  scope  of  this  study. 
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Figure  24  CORRELATION  OF  DATA  WITH  WEBER  NUMBER 


3.4  STRIPPING  MODE 

An  analysis  for  the  mass  loss  during  stripping  has  been  given  by  Taylor.**  Taylor 
observed  that  the  motion  of  the  air  over  the  surface  of  the  drop  tends  to  drag 
fluid  along  with  it,  and  he  assumed  that  such  fluid  might  be  spilled  over,  or 
separated,  into  the  wake  of  the  flow  somewhere  near  the  equator  of  the  drop, 
where  lta  velocity  is  greatest.  His  analysis  did  not  go  into  the  details  of 
the  spillage  process,  but  simply  required  that  the  accelerated  fluid  (or  some 
fraction  thereof)  in  the  boundary  layer  on  the  surface  of  the  drop  be  ripped  off 
in  some  manner  when  it  reaches  the  equator.  The  problem  then  reduced  to  calcu¬ 
lating  the  mass  flow  in  the  boundary  layer,  which  varies  as  a  function  of  the  drop 
diameter,  speed,  etc. 

There  la  an  error  in  the  Taylor  analysis  involving  the  position  of  the  parameter  A 
in  the  final  result.  Ranger  and  Nlcholls*  have  corrected  Taylor's  error,  and 
they  report  adequate  agreement  of  the  theory  with  their  data.  Although  the 
deceleration  of  the  drop  can  be  included  in  the  analysis,  we  shall  start  by 
applying  Ranger  and  Nlcholls'  theoretical  result  to  the  case  of  no  deceleration. 
The  stripping  analysis  then  predicts  the  following  law  of  variation  of  the  mass 
loss  rate: 

TT"  ‘  ("l»5/a  (j  \  ^  «2a)1/2  (28) 
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where 


inn  x1  5 

V'l  n ) 
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For  a  spherical  drop,  the  mass  is  related  to  the  drop  diameter  by  the  equation 


T  nD- 


(30) 


We  then  have,  in  terms  of  the  nondimensional  time, 
/  m  \1/2  d  /  m  \  f  n  A 

\  mo  /  dT  \  mo  /  V  ”  n  u2  Do  ’ 

which  has  the  solution 


(31) 
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(32) 


Hence,  the  relevant  parameter  for  the  stripping  mode  would  be  the  "mixed  Reynolds 
number" 


I' 2  u  2  Do 
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The  Mach  3  data  shown  in  Figure  17  present  the  conditions  under  which  stripping 
apparently  progressed  farthest  until  the  onset  of  large  instabilities.  An  extra¬ 
polation  of  these  data  suggests  that,  had  stripping  continued,  the  drop  mass 
would  have  vanished  at  between  180  and  200  msec,  i.e.,  at  a  value  of  T  of  about 
3.5.  The  value  of  R  for  these  conditions  is  approximately  94,000.  If  we  sub¬ 
stitute  this  value  into  Equation  (33)  and  compute  the  nondimensional  time 
corresponding  to  m  ■  0,  we  obtain  T  ■  47,  more  than  an  order  of  magnitude 
greater  than  the  experimental  breakup  time.  It  Beems  unlikely  that  any  simple 
correction  to  the  theoretical  result  would  be  sufficient  to  overcome  such  a 
large  discrepancy.  We  are  therefore  led  to  the  conclusion  that  the  analysis 
developed  by  Taylor  and  corrected  by  Ranger  and  Nlcholls  does  not  represent  a 
valid  theoretical  interpretation  of  the  observed  stripping  of  mass  from  water 
drops  at  high  relative  velocities;  the  mechanism  of  stripping  mist  be  roughly  an 
order  of  magnitude  more  effective  than  that  suggested  by  these  authors. 

We  shall  not  attempt  to  provide  a  new  analytical  representation  of  the  stripping 
mode  in  this  report,  but  rather  concentrate  on  obtaining  a  simplified  correlation 
of  the  experimental  data  which  will  be  adequate  for  flight  predictions. 


-33 


Further  consideration  of  the  Mach  3  data,  as  well  as  the  Mach  6  results  (Fig¬ 
ure  18),  indicates  that  stripping  initially  reduces  the  mass  at  a  fairly  slow 
rate  and  subsequently  becomes  more  effective;  i.e.,  dm/de  is  small  at  early  time. 
At  the  other  extreme,  i.e.,  when  stripping  is  well  established,  we  may  postulate 
that  the  mass  removal  rate  is  proportional  to  the  drop  mass  remaining,  which 
would  suggest  that  dm/dt  *  0  att  -  t  ,  where  the  subscript  s  indicates  breakup 
time  due  to  stripping.  Hence  the  curve  m(t )  may  be  expected  to  have  a  small 
slope  both  at  small  time  and  near  ts  .  A  simple  curve  satisfying  these  require¬ 
ments  is 
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(34) 


Using  a  value  of  ts  -  3.5,  this  curve  is  shown  along  with  the  experimental  mass 
histories  in  Figures  25,  26,  and  27.  The  agreement  with  the  data  is  good  at  all 
three  test  conditions.  Note  that  no  theoretical  claim  can  be  made  concerning 
Equation  (34).  It  is,  however,  a  physically  reasonable  and  accurate  representa¬ 
tion  of  the  experimental  data.  Note  also  the  apparent  constant  value  of  Tb  . 

We  have  plotted  the  lower  Weber  number  data  of  Engel, 2  Nicholson,3  and  Ranger  and 
Nicholls4  versus  R  in  Figure  28.  We  have  also  included  our  Mach  3  print  -ifth  its 
estimated  Tb  due  to  stripping  of  3.5.  Although  there  appears  to  be  a  slight  posi¬ 
tive  dependence  of  Tb  upon  R,  the  large  scatter  makes  this  questionable.  In  fact. 
Nicholson  interpreted  his  own  data  as  Indicating  a  constant  value  fir  rb  of  3.5. 
The  stripping  analysis  of  Equation  (32)  is  therefore  open  to  serious  question, 
since  no  dependence  on  the  mixed  Reynolds  number  can  be  inferred  fiom  the  data. 

With  regard  to  choosing  a  working  value  for  Tb  due  to  stripping,  we  are  inclined 
tentatively  to  conclude  that  it  is  3.5  and  constant.  This  value  lies  between  the 
results  of  Engel  and  Ranger  and  Nlcholls,  is  quoted  by  Nicholson  on  the  basis  of 
his  more  quantitative  definition  of  drop  breakup,  and  is  compatible  with  our 
own  estimated  value  on  the  basis  of  the  X-ray  data. 


Ftywe  25  MOP  MASS  HISTORIES  AT  MACH  S 
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4.0  DROP  BREAKUP  CALCULATIONS  UNDER  FLIGHT  CONDITIONS 


Employing  the  experimental  and  theoretical  criteria  for  drop  breakup  reported 
above,  representative  prelimina,  calculations  have  been  made  to  assess  the 
importance  of  drop  breakup  under  realistic  reentry  conditions.  The  results 
suggest  that  the  primary  mode  in  the  stagnation  region  is  catastrophic  breakup, 
while  the  stripping  mode  is  dominant  in  the  downstream  region  of  the  vehicle 
shock  layer. 

Drop  breakup  calculations  are  naturally  divided  into  three  parts.  In  the  first 
part,  the  experimental  correlation  of  drop  deceleration  in  the  shock  tube  is 
expressed  in  vehicle-fixed  coordinates,  and  the  nondimens ional  impact  time  is 
then  computed  as  a  function  of  drop  diameter.  In  the  second  part,  the  drop 
breakup  criteria  are  used  to  ascertain  the  appropriate  nondimensional  breakup  time 
as  a  function  of  drop  diameter.  Finally,  the  two  times  are  compared  in  order  to 
determine  the  fraction  of  drop  mass  entering  the  shock  layer  which  does  not  reach 
the  vehicle  surface  as  a  function  of  drop  diameter,  and  this  function  is  integrated 
over  the  rainstorm  drop  mass  distribution  function  to  yield  the  net  mass  fraction. 

Two  cases  have  been  considered:  the  stagnation  region  of  a  blunt-nosed  body 
and  the  conical  region  far  downstream  of  the  nose.  In  the  stagnation  region,  a 
reasonable  estimate  is  obtained  by  considering  only  the  breakup  of  drops  which 
traverse  the  stagnation  streamline.  Although  drops  which  traverse  other  parts 
of  the  shock  layer  will  experience  somewhat  different  flow  conditions,  the 
stagnation  streamline  provides  representative  conditions  and  is  simple  to  analyze. 
Furthermore,  a  direct  analogy  relates  the  stagnation  streamline  problem  to  the 
experimental  conditions  in  the  shock  tube:  In  a  high-speed  shock  layer,  the 
flow  conditions  along  the  stagnation  streamline  are  virtually  uniform  and  hence 
similar  to  the  uniform-flow  region  behind  the  shock  wave  in  the  shock  tube.  In 
the  downstream  region,  the  disturbance  caused  by  the  blunt  nose  tends  to  die  out, 
and  the  flow  field  becomes  approximately  conical.  High-speed  conical  shock  layers 
are  also  characterized  by  approximate  uniformity,  and  similarity  with  experimental 
conditions  can  be  established  if  the  relative  drop  motion  in  the  direction  per¬ 
pendicular  to  the  vehicle  surface  is  considered. ^  Although  it  is  possible  to 
consider  more  complicated,  non-uniform  flow  fields,  this  is  beyond  the  scope  of 
the  present  effort. 

In  shock-tube  coordinates,  the  distance  between  the  shock  wave  and  the  drop  is 
given  by 

s  *  “s  1  "  *  >  (35) 

where  x  is  the  distance  the  drop  has  moved  from  its  original  position.  In 
nondimensional  variables,  the  experimental  drop  acceleration  correlation  is 

X  -  0.8  T2  .  (36) 

Since 
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Equation  (35)  can  be  written 
S  -  B  T  -  0.8  T2  , 

where 


b„r  pJp  L- 


This  equation  has  the  solution 


(38) 


(39) 


(40) 


In  relating  this  equation  to  the  case  of  flight  through  the  atmosphere,  we 
interpret  us  as  the  flow  velocity  component  normal  to  the  shock  wave:  in  the 
stagnation  region,  the  free-stream  velocity  ;  in  the  conical-flow  region,  the 
velocity  sin  0  ,  where  0s  is  the  inclination  of  the  shock  wave.  In  what 

follows,  we  shall  be  concerned  only  with  conditions  at  impact  on  the  vehicle 
surface.  In  the  stagnation  region,  the  nondimens ional  distance  S  is  then 
interpreted  as  A/DQ  ,  where  A  is  the  shock  standoff  distance,  related  to  the 
nose  radius  R^  and  the  shock  density  ratio  t  by  the  high-speed  correlation^ 

A  ~  0.78  (  Rn.  (41) 

In  the  conical-flow  region,  S  is  interpreted  as  the  perpendicular  distance  from 
the  shock  wave  to  the  body  L  .  »  Rj  sin  ( 0S  -  0C ) ,  where  Rj  is  the  distance  from  the 
cone  vertex  at  which  the  drop  strikes  the  vehicle  surface.  We  can  evaluate  L  by 
the  high-speed  correlation^ 

L  ?  0.6  r  es  R; .  (42) 


Here  0Q  is  the  inclination  of  the  vehicle  surface. 
(40)  then  gives  the  impact  time  in  the  form 
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For  the  two  regions,  Equation 


(43) 


where 

ZN  »  2.5  (1  -<)2b^ 

Zc  -  l.92(l-02b„»s.  (44) 
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It  is  therefore  seen  that,  for  given  flight  conditions  and  for  a  given  ,  the 
nondimens ional  impact  times  are  functions  only  of  the  ratios  RN/DQ  and  Rj  / D0  . 


We  have  considered  three  altitudes  (0,  15,  and  30  kft)  and  three  flight  Mach 
numbers  (5,  10,  and  20).  Figures  29,  30,  and  31  show  the  impact  time  in  the 
nose  region  as  a  function  of  RN/DQ  f°r  these  conditions.  For  impact  downstream 
on  the  vehicle,  we  have  considered  three  cone  angles  (  0Q  ■  5,  10,  and  15  degrees). 
Instead  of  plotting  the  corresponding  impact  time  from  Equation  (43)  for  all  27 
cone  conditions,  we  note  that  for  values  of  B  greater  than  about  30,  the  accelera¬ 
tion  term  0.8  Tj  is  very  small  and  can  be  dropped  from  Equation  (38);  the  impact 
time  is  then  simply  given  by  the  relationship 
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where 


C  ,  0.6 


(46) 


For  those  cases  in  which  this  relationship  is  valid,  the  factor  C  is  tabulated 
in  the  following  table. 
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0.525 

10 
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1.71 
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1.05 

10 

1.38 

1.07 

15 

5 

2.72 

2.16 

1.66 

The  experimental  breakup  criteria  which  we  have  used  are  as  follows:  If  the 
impact  t*me  T;  is  greater  than  the  breakup  time  Tb ,  the  drop  mass  at  impact  is 
taken  to  be  zero.  If  the  impact  time  is  less  than  the  breakup  time,  the  drop 
mass  at  impact  is  given  by  the  stripping-mode  correlation 
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Figure  31  TIME  TO  IMPACT  STAGNATION  POINT ,  h  =  30  kft 


where  n^,  is  the  initial  mass.  The  breakup  time  is  taken  to  be  the  smaller  of  3.5 
(stripping  mode)  and  45 W_1/^(catastrophic  mode),  where 
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'V  sin  2  6  , 


(48) 


It  should  be  noted  that  the  range  of  experimental  data  which  led  to  the  latter 
correlation  for  catastrophic  breakup  is  exceeded  somewhat  under  flight  conditions 
in  the  stagnation  region. 

We  have  only  considered  the  case  Mm  =*  10  at  15  kft  altitude  in  the  detailed  breakup 
calculation  below.  At  that  altitude,  there  is  a  peak  in  the  rainstorm  water  con¬ 
tent  for  the  particular  model  which  we  chose  (described  below),  so  that  these 
conditions  are  representative  of  a  severe  rain  environment.  We  chose  representa¬ 
tive  nose  radii  of  0.02,  0.1  and  0.5  foot,  and  in  the  conical-flow  region  we  chose 
a  10-degree  half-angle  cone  with  Rj  *  1  and  5  feet.  (These  values  essentially 
designate  the  axial  location  of  the  impact  points.)  The  standoff  distance  \  for 
these  conditions  is  about  0.002,  0.01,  and  0.05  foot,  respectively  (A/RN*>  0.0975, 

<  =  0.125). 

Figure  32  shows  the  impact  times  for  the  designated  conditions  as  functions  of 
the  dimensionless  distances.  Figures  33  and  34  show  the  impact  and  breakup 
times  as  functions  of  the  initial  drop  diameter  (  D0  )  in  millimeters  for  the  four 
cases.  Figure  33  shows  that  the  breakup  time  in  the  nose  region  is  essentially 
dominated  by  the  catastrophic  mode.  We  observe  that  in  the  nose  region,  drops 
with  diameters  less  than  0.17  mm  do  not  reach  the  surface  of  the  0.1-foot-radius 
nose,  while  for  the  0.5-foot-radius  nose  the  breakup  effect  is  much  more  important: 
drops  smaller  than  1.4  mm  diameter  are  dissipated. 

In  the  conical-flow  region  (see  Figure  34),  although  the  catastrophic  mode  affects 
the  results  for  the  5-foot  station  slightly,  the  stripping  mode  largely  controls 
the  amount  of  mass  reaching  the  surface.  At  the  1-foot  station,  drops  with 
diameters  smaller  than  0.12  mm  are  completely  dissipated  before  they  reach  the 
surface,  while  at  the  5-foot  station  the  corresponding  diameter  is  0.7  mm.  We 
chose  as  a  rainstorm  model  a  typical  summer  storm  referred  to  in  the  Air  Force 
Handbook  of  GeophyslcslO  as  Case  1.  This  case  is  believed  to  be  representative 
of  a  summer  rain  in  temperate  latitudes,  and  to  come  close  to  simulating  widespread, 
tropical  rainB.  It  was  derived  by  combining  empirical  and  theoretical  data  in¬ 
dicating  approximate  weather  conditions  in  this  type  of  Btorm.  The  maximum 
updraft  was  taken  to  be  0.2  m/sec  yielding  a  precipitation  rate  of  0.2  in. /hr. 

Such  a  rate  would  be  exceeded  in  Washington,  D.C.,  about  40  hours  in  an  average 
year  and  could  probably  be  referred  to  as  a  medium  rain.  The  precipitation  rate 
will  not  enter  directly  into  these  calculations,  however;  the  results  will  be 
given  in  terms  of  the  mass  fraction  of  the  incoming  rain  that  reaches  the  vehicle 
surface,  so  that  the  impinging  mass  rate  per  surface  area  of  the  vehicle  must  be 
determined  by  multiplying  by  the  rain  concentration,  the  vehicle  velocity,  and 
the  surface  inclination.  At  15  kft  altitude,  the  rain  concentration  for  this 
storm  is  about  1  gm/m^,  giving  an  incoming  mass  rate  of  about  0.02  slug/ft2-sec. 
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NONDIMENSIONAL  IMPACT  TIME 
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Figure  32  NOSE  AND  CONE  IMPACT  TIMES 
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Figure  33  IMPACT  AND  BREAKUP  TIMES  FOR  NOSE 
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Figure  34  IMPACT  AND  BREAKUP  TIMES  FOR  CONE 


Reference  10  (page  6-9)  gives  an  expression  for  the  raindrop  sise  distribution 
function  which  can  be  converted  into  the  following  ness  distribution  function: 


Here  pdD  is  the  ness  fraction  of  drops  with  disasters  between  D  and  D  +  dD ,  and 
the  reference  diameter  Op  is  related  to  the  reference  diameter  "  D.  "  in  Reference 
10  by 


Dr  -  "D0"/3.67 


(50) 


This  distribution  function  is  shown  plotted  in  Figure  35 %  The  reference  diameter 
Dr  is  a  weak  function  of  the  precipitation  ratelu;  in  this  case,  DR  -  0.255  mm. 


It  should  be  noted  that  this  model  storm  takes  into  account  only  precipitating 
rain  and  neglects  the  cloud  water  content,  which  in  general  is  significantly 
greater  than  the  precipitating  water  content.  The  cloud  water  is  also  distributed 
over  much  smaller  particles,  so  that  the  shock  layer  effect  on  the  cloud  water 
should  be  more  pronounced  than  on  the  precipitating  water. 
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The  mass  fraction  of  drops  with  diameters  between  0  and  a  given  diameter  O^can  be 
determined  by  integrating  Equation  (49)  between  the  respective  limits: 


fxdD  =  P  (4,  Dj)  , 


(51) 


where 

P  (4,  Dj) 


1 

D,  +  - 
1  2 


-  2 
(D  j) 


-D 

e 


1  , 


Di  =  Dj/Dr 


(52) 


We  have  computed  the  mass  fraction,  Mj  , 
of  drop  mass  which  strikes  the  surface. 

In  the  nose  region  this  can  be  accom¬ 
plished  with  reasonable  accuracy  by  eval¬ 
uating  Equation  (51)  in  terms  of  the 
diameter  Dj  of  the  smallest  particle  to 
reach  the  surface,  previously  cited;  then 
Mj  -  1  -  Mj  .  The  impact  times  are  gen¬ 
erally  so  low  that  there  is  no  signifi¬ 
cant  effect  of  stripping  on  the  remaining 
drops.  In  the  conical-flow  region,  how¬ 
ever,  such  an  estimate  would  disregard 
the  fact  that  a  significant  fraction  of 
the  drops  which  do  reach  the  surface  have 
had  part  of  their  mass  stripped  off. 

The  corresponding  mass  fraction  lost  must 
be  computed  by  applying  Equation  (47)  to 
the  mass  distribution  function  and  inte¬ 
grating  over  the  remaining  drops.  The 
total  impinging  mass  fraction  then  becomes 


(53) 


where  D;  is  the  drop  diameter  for  which  Tj -  Tb ,  TO  facilitate  this  computation, 
a  digital  computer  program  has  been  written. 


-44- 


The  results  of  the  mass  removal  calculations  are  summarized  in  the  following  table 
for  both  the  nose  and  cone  regions. 


10,  h  -  15  kft,  0Q  10  degrees 


Critical  Drop  Diameter, 
D,  (mm) 

Impinging  Mass 
Fraction, M; 

Nose 

Rn  =  0.02  ft 

0.02 

1.00 

rn  =  0.1  ft 

0.17 

0.97 

RN  *  0.5  ft 

1.4 

0.22 

Cone 

*  1  ft 

0.12 

0.92 

Ri  =  5  ft 

0.75 

0.29 

Clearly  for  the  0.02-  and  0.1-foot-radius  nose  and  1-foot  station  on  the  cone, 
the  drops  lose  virtually  none  of  their  mass,  and  the  impinging  mass  on  the  surface 
at  these  points  will  comprise  nearly  all  the  mass  which  has  entered  the  shock 
layer  at  the  corresponding  points.  For  the  0.5-foot-radius  nose,  very  little 
mass  strikes  the  vehicle  surface  at  the  stagnation  point,  and  the  vehicle  shock 
layer  is  78  percent  efficient  in  protecting  the  surface  from  the  rain.  At  the 
5-foot  station  on  the  cone,  71  percent  of  the  entering  rain  mass  will  be  dissipated 
through  the  mechanism  of  drop  breakup. 

Figure  36  shows  details  of  the  impinging  drop  mass  fraction  as  a  function  of  drop 
diameter  (in  microns)  for  the  0.02-foot-radius  nose.  It  may  be  noted  that  for 
this  small  nose  radius,  the  stripping  mode  becomes  important  for  drop  diameters 
very  near  the  cutoff  diameter — 20  microns — for  catastrophic  breakup.  Since  our 
predictions  are  based  on  correlations  of  data  obtained  for  much  larger  drop  sizes, 
additional  experimental  work  would  be  required  to  verify  this  interesting  conclu¬ 
sion  and  to  pin  down  the  critical  drop  breakup  parameters  for  very  small  droplets. 

Figures  37  through  39  show  drop  Impact  velocities  and  angles  as  functions  of  drop 
diameter,  D0  .  The  drop  impact  velocity  components  normal  to  the  body  surface 
were  obtained  by  differentiating  the  experimental  drop  acceleration  correlation. 
Equation  (36),  The  impact  angles  were  obtained  from  the  impact  velocity  component 
normal  to  the  surface  together  with  the  component  tangential  to  the  surface;  the 
latter  was  assumed  unchanged  during  passage  through  the  shock  layer. 

It  can  be  concluded  that  drop  breakup  effects  are  potentially  important  in 
high-speed  shock  layers.  Depending  on  the  flight  conditions  and  vehicle  con¬ 
figuration,  a  large  fraction  of  the  entering  drop  mass  may  be  dissipated  in  the 
shock  layer,  never  reaching  the  vehicle  surface.  In  the  stagnation  region,  the 
catastrophic  breakup  mode  la  extremely  important.  Correlations  of  the  experimental 
data  presented  in  this  report  can  be  used  to  make  detailed  predictions  of  raindrop 
breakup  effects  for  reentry  vehicles. 
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Figure  36  IMPINGING  DROP  MASS  FRACTION  VERSUS  DIAMETER  FOR 
0.02-FOOT-RADIUS  NOSE 


Figure  37  IMPACT  VELOCITY  FOR  NOSE 
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Figure  38  NORMAL  VELOCITY  AT  IMPACT  FOR  CONE 


Figure  39  IMPACT  ANCLE  FOR  CONE 
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5.0  SUMMARY  AND  CONCLUSIONS 


An  oxperimental  program  has  been  carried  out  to  determine  the  modes  of  raindrop 
breakup  in  air  after  passage  of  a  strong  shock  wave.  Relevant  dimensionless 
parameter 3  have  been  determined,  and  the  experimental  data  have  been  used  to 
determine  raindrop  breakup  in  the  shock  layer  of  a  reentering  vehicle  passing 
through  a  rainstorm. 

The  raw  experimental  data  consisted  of  both  shadowgraphs  and  X-ray  photographs  of 
disintegrating  waterdrops  in  a  shock  tube.  Data  were  taken  for  a  range  of  drop 
diameters  from  0.5  to  2.5  mm,  for  ambient  pressure  levels  from  140  to  760  torr, 
and  for  shock  Mach  numbers  from  3  to  11.  The  X-ray  photographs  were  taken  by 
doping  the  drops  to  saturation  with  lead  acetate.  The  motivation  behind  using 
X-ray  photography  was  to  obtain  reliable  data  on  the  drop  mass  remaining  at  late 
stages  in  the  disintegration  process,  when  the  cloud  of  mist  previously  torn 
from  the  drop  totally  obscures  rht  drop  contour  on  a  shadowgraph.  In  the  regime 
dominated  by  catastrophic  breakup,  the  X-ray  results  show  qualitative  agreement 
with  conventional  spark  shadowgraph  data,  but  indicate  about  half  of  the  times 
to  breakup  obtained  from  the  shadowgraph  method.  Determination  of  the  mass 
remaining  in  the  drops  was  accomplished  by  processing  the  X-ray  negatives  with 
a  microdensitometer,  and  then  integrating  the  resulting  output  digitally. 

Examination  of  the  shadowgraphs  demonstrated  conclusively  that  material  was  being 
stripped  from  the  drops  continuously  at  all  test  conditions.  At  the  more  severe 
test  conditions,  a  second  "catastrophic"  mode  of  breakup  became  dominant,  and 
the  drop  was  shattered  before  all  its  material  could  be  removed  by  stripping. 

The  catastrophic  mode  of  breakup  is  related  to  the  unstable  growth  of  waves  on 
the  front  surface  of  the  drop  caused  by  the  rapid  acceleration  of  the  drop  by  the 
air  stream.  In  this  mode,  surface  waves  are  observed  to  grow  rapidly  until  their 
amplitude  is  comparable  to  the  dimensions  of  the  drop,  at  which  point  the  drop 
is  torn  apart. 

An  additional  effect  which  is  superimposed  on  these  two  modes  of  breakup  is  the 
initial  flattening  of  the  drop  caused  by  the  imposed  pressure  differential  between 
the  poles  and  the  equator.  Because  the  deformation  of  the  drop  may  affect  the 
breakup  process,  shadowgraphs  were  taken  at  the  early  stages  of  the  process  to 
determine  the  rate  of  deformation.  The  deformation  of  the  drop  also  affects 
its  trajectory,  since  the  dreg  coefficient  cannot  be  considered  to  remain 
constant.  In  making  flight  predictions,  it  was  decided  to  rely  on  the  experimental 
data  for  Information  on  the  drop  trajectory,  and  these  date  were  therefore  reduced 
and  correlated. 

Earlier  investigators  hava  found  that  the  data  on  drop  trajectories  could  be 
correlated  by  the  relationship 

X  0.8  T2 

where  X  and  T  are -  dimensionless  shock-tube  distance  and  time,  respectively.  The 
data  on  initial  drop  datamation  wars  at  varianca  with  Burgers'  analytical  result, 
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but  could  be  correlated  by  a  simplified  analytical  model  Involving  the  assumption 
of  a  disk-shaped  drop.  The  failure  of  the  early-time  analysis  to  predict  the 
drop  deformation  is  apparently  due  to  the  fact  that  the  range  of  applicability  of 
the  early-time  assumption  is  vanishingly  small. 

The  catastrophic-mode  data,  obtained  from  the  X-ray  photographs,  could  be 
correlated  by  the  formula 

Tfa  45  W~  1/4 


where  Tb  is  the  time  to  breakup  and  W  is  the  Weber  number.  This  variation  is  in 
accord  with  the  small-disturbance  theory  for  unstable  waves. 


It  was  decided  to  take  the  dimensionless  time  to  breakup  for  the  stripping  mode 
equal  to  3.5,  based  on  the  experimental  data.  This  was  done  with  some  reluctance, 
since  the  theory  for  stripping  predicts  a  variation  for  Tb  as  the  1/2-power  of 
the  mixed  Reynolds  number  defined  by  Equation  (31).  There  was  however  insufficient 
data  to  justify  such  a  conclusion,  and  the  available  experimental  data  are 
reasonably  correlated  by  the  value  3.5.  The  data  for  the  mass  loss  during  the 
stripping  process  could  be  correlated  by  the  formula 


where  m0  is  the  initial  mass  of  the  drop  and  m  is  the  mass  of  the  drop  at  the 
dimensionless  time  T . 

Employing  these  correlations,  calculations  of  the  dissipation  of  raindrop  mass  in 
the  shock  layer  of  a  reentering  vehicle  have  been  carried  out.  The  relationship 
between  the  experimental  setup  and  the  passage  of  drops  through  the  blunt-nose 
region  and  through  the  downstream  conical-flow  region  has  been  determined,  and 
the  appropriate  transformations  have  been  made.  Dimensionless  impact  times  have 
been  calculated  for  a  number  of  flight  conditions.  The  specific  case  of  drop 
breakup  during  flight  through  the  15  kft  level  of  a  representative  rainstorm  has 
been  considered.  It  was  determined  that  for  a  0.1-foot  nose  radius  and  for  the 
1-foot  station  on  a  10-degree  cone,  essentially  no  mass  is  lost  by  the  raindrops 
before  they  strike  the  surface  of  the  vehicle.  For  a  0.5-foot  nose  radius,  most 
of  the  incoming  drop  mass  is  dissipated,  and  the  dominant  mechanism  for  breakup 
is  the  catastrophic  mode.  At  the  5-foot  station  on  the  10-degree  cone,  about 
half  of  the  incoming  drop  mass  is  dissipated;  the  dominant  mechanism  here  is  the 
stripping  mode. 

We  conclude  that  computations  of  raindrop  breakup  in  vehicle  shock  layers  can  be 
carried  out  under  realistic  conditions  using  the  data,  techniques,  and  correlations 
developed  in  this  report.  Two  modes  of  drop  breakup  have  been  found  to  be  import¬ 
ant  for  high-speed  flow.  In  the  blunt-nose  region  the  catastrophic  node  of 
breakup  appears  to  be  dominant,  whils  in  the  downstream  flow  region  the  stripping 
mode  is  the  primary  mechanism. 
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APPENDIX 


TABLE  OF  EXPERIMENTAL  DATA 


Shadowgraph  Data:  Test  liquid  was  distilled  water  except  where  asterisk 
indicates  saturated  lead  acetate  solution. 


Nominal 

Initial 

Drop 

Diameter 

(mm) 

Nominal 

Shock 

Mach 

Number 

Initial 

Pressure 

(torr) 

Normalized 

Drop 

Diameter 

Dimensionless 
Time  after 
Shock  Passage 

1.0 

1.00 

0 

1.0 

Hi 

■ 

1.035 

1.0 

■1 

350 

1.12 

1.0 

350 

1.21 

■ 

1.0 

140 

1.30 

0.33 

140 

1.62 

0.46 

6 

350 

1.37 

0.29 

6 

350 

1.18 

0.14 

0.5 

9 

350 

1.52 

0.34 

0.5 

1  1 

140 

2.84 

1.43 

0.5 

EH 

140 

3.20 

1.34 

0.5 

140 

3.50 

1.88 

0.5 

i 

140 

3.25 

2.50 

0.5 

i 

140 

Gone 

2.75 

0.5 

H 

140 

1.12 

0.15 

2.5 

WSm 

350 

1.47 

0.33 

2.5 

9 

350 

1.18 

0.22 

2.5 

12 

90 

1.54 

0.41 

2.5 

3 

350 

1.23 

0.20 

2.5 

3 

350 

1.39 

0.34 

2.5 

3 

350 

1.27 

0.27 

1* 

6 

350 

1.63 

0.40 

1* 

6 

350 

3.57 

1.28 

1* 

6 

350 

4.01 

3.25 

1* 

6 

350 

3.40 

2.56 

1* 

6 

350 

2.38 

4.02 

1* 

6 

350 

4.14 

3.68 

1* 

6 

350 

2.49 

0.78 

1* 

6 

350 

1.44 

0.30 

1* 

6 

350 

2.99 

2.31 
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X-Ray  Data:  Test  liquid  was  saturated  lead  acetate  solution 
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Initial 

Drop 

Diameter 

(mm) 

Nominal 

Shock 

Mach 

Number 

Initial 

Pressure 

(torr) 

Normalized 

Drop 

Diameter 

Dimensionless 
Time  after 
Shock  Passage 

1 

6 

mm 

0.99,1.00 

0 

1 

6 

WMM 

0.62,0.64 

1.38 

2.12 

6 

350 

0.45 

1.95 

2.12 

6 

350 

0.23 

1.99 

2.12 

6 

350 

0.99,0.89, 

0.96 

0.31 
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2.12 

3 
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0.92 

0.76 

2.12 

3 

760 

0.81 

0.45 

2.12 

3 

760 

0.48 

1.28 

2.12 

3 

760 

0.66,0.74 

1.46 

2.12 

3 

760 

0.34 

2.58 

2.12 

3 

760 

0.82 

1.59 

2.12 

3 

760 

0.37 

1.68 

2.12 

3 

760 

0.65 

1.19 

2.12 

11 

141 

1.01,0.94 

0.15 

2.12 

11 

141 

1.05 

0 

2.12 

11 

141 

0.92 

1.32 

2.12 

11 

141 

0.89 

1.19 

2.12 

11 

141 

0.76 

1.07 

2.12 

11 

141 

0.77 

0.57 

2.12 

11 

141 

0.71,0.80 

0.94 
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